Evaluation of Finite-Rate Gas/Surface Interaction 


Models for a Carbon Based Ablator 


Yih-Kanq Chen* 

NASA Ames Research Center, Moffett Field, CA 94035-1000 

and 

Tahir GdkceiV 

ERC, Inc., NASA Ames Research Center, Moffett Field, CA 94035-1000 

Two sets of finite-rate gas/surface interaction model between air and the carbon surface are studied. The 
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surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature 
diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface 
recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The 
material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and 
Ablation Program, which predicts charring material thermal response and shape change on hypersonic 
space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line 
Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc- 
jet facilities with heat fluxes ranging from 45 to 1100 w/cm 2 are performed and compared with data for 
model validation. The ablating material used in these arc -jet tests is Phenolic Impregnated Carbon Ablator. 
Computational predictions of surface recession and shape change are in good agreement with measurement 
for arc -jet conditions. 
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= normalized ablation rate 
= mass fraction 
= diffusion coefficient, m 2 /s 
= enthalpy, J/kg 

= forward reaction rate, (units vary) 

= thermal conductivity of translational temperature, W/m-K 
= thermal conductivity of vibrational temperature, W/m-K 
= molecular weight, kg/kmol 
= mass flux, kg/m 2 -s 
= pressure. Pa 

= in-depth conductive heat flux at surface, W/m 2 
= convective heat flux at surface, W/m 2 
= radiative heat flux at surface, W/m 2 
= corner radius of the model, cm 
= nose radius of the model, cm 
= universal gas constant, 8.3143 J/mol-K 
= temperature, K 

= maximum allowed change on surface temperature, K 


T t = translational temperature, K 


T v = vibrational temperature, K 

= environment temperature, K 

v = mass injection velocity, m/s 

V = thermal velocity, m/s 

w = source term of gas-surface interactions, kmol/m 2 -s 
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a 


= absorptance 


£ 

= emissivity 

P 

= total density, kg/m 3 

0 

= active site density at the surface, mol/m 2 

cr 

= Stefan-Boltzmann constant, W/(m 2 -K 4 ) 

subscripts 

c 

= char 

8 

= pyrolysis gas 


i = gas species 

w = wall 


I. Introduction 

The function of an atmospheric entry vehicle’s Thermal Protection System (TPS) is to protect underlying materials 
from the intense aerothermal heating associated with hypersonic entry. Typical ablating TPS materials pyrolyze and 
ablate at high temperature for mass-efficient rejection of the aerothermal heat load. For design and sizing of ablating 
TPS materials, it is imperative to have reliable numerical procedures which can accurately compute the vehicle 
aerothermal environment, the material surface ablation, and the material internal thermal response. It has been 
demonstrated that accurate prediction of ablative heat flux requires a fluid-solid shape change coupling 
simulation. [1,2] In our early fluid-solid shape change coupling computations, the ablative carbon surface was assumed 
to be at chemical equilibrium. Chemical equilibrium is a good assumption for many space entry applications, but it 
may not be valid for all conditions. Furthermore, the aerothermal environments computed by Data Parallel Line 
Relaxation method (DPLR)[3] were for a fully catalytic non-ablating surface. Thus, an engineering correlation with a 
blowing reduction parameter had to be introduced in the Two-dimensional Implicit Thermal -response and AblatioN 
program (TITAN)[4] simulation to take into account the effect of mass injection on reduction of convective heat flux. 
A mass transfer coefficient also had to be defined based on the heat transfer coefficient under the assumption of a 
relatively thin boundary layer for the computation of char recession rate. In our recent studies, integrated fluid-material 
response analyses with finite-rate gas-surface interactions were performed. The coupling procedure, considering 
surface thermal chemistry and shape change, was achieved by using either explicit or implicit scheme [5,6], However, 
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for explicit coupling simulations with surface thermal chemistry included, numerical oscillation of predicted surface 
temperature was observed at conditions with relatively high surface heat fluxes. The time step size for updating the 
interface boundary conditions had to be reduced to minimize the numerical oscillation. Consequently, the total 
computational time required to complete a simulation increased significantly. Additionally, for some explicit coupling 
simulations, the numerical oscillation of surface temperature became so severe that a converged solution was not 
obtained. An implicit approach was proved to be more efficient for performing fluid-solid surface thermal chemistry 
coupling for the high surface heating arc -jet environments and the trajectory-based simulations, such as that of Stardust 
Earth Reentry Capsule. 

The air-carbon surface interactions studied in our previous work were mostly in the so called diffusion -controlled 
oxidation or low sublimation regime. For surface temperatures below the diffusion-controlled regime, the net carbon 
oxidation rates have to be determined by both the forward and backward reaction rates. This is the so called rate -(or 
kinetic) controlled carbon oxidation. The engineering model developed by Park’s[7] one-way surface chemical 
reaction may be inadequate to precisely predict the carbon recession if the backward surface reaction rates are 
significant. The two-way surface reaction model developed by Zhluktov and Abe [8] were shown to be far under- 
predicting both the carbon recession rate and surface heat flux even if the surface temperature is in the diffusion 
controlled regime. [5] 

The primary objective of this work is to enhance the existing finite -rate carbon/air surface chemistry models to 
improve the prediction accuracy of surface recession and surface heat flux for carbon based ablators. These revised 
gas/surface interaction (GSI) models will be applicable to carbon oxidation at both rate and diffusion controlled 
conditions as well as carbon sublimation at a much elevated surface temperature. The DPLR/TITAN coupled 
simulation system is used in the present work to predict the ablator surface recession and shape change. The boundary 
conditions with general non-equilibrium finite-rate chemistry for gas/surface interactions implemented in the DPLR 
code by MacLean[9] is adopted to predict char mass injection rate. The surface species mass balance is thus performed 
in DPLR, and the surface energy balance is performed in TITAN. The hot-wall ablating convective heat flux is directly 
computed in DPLR based on the surface temperature and pyrolysis gas injection rate computed in TITAN. Coupled 
fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc -jet facilities with 
heat fluxes ranging from 45 to 1100 w/cm 2 are performed using the GSI models and their rates proposed in this work. 
The ablating material used in these arc -jet tests is a Phenolic Impregnated Carbon Ablator (PICA). [10] Computational 
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predictions of surface recession, and shape change for arc -jet test conditions of Small Probe Reentry Investigation for 
TPS Engineering (SPRITE) [12] test article are also performed and compared with the measurement. 

II. Gas/Surface Interactions 

Two sets of non-equilibrium finite-rate gas-surface interaction models between air and carbon surface are proposed 
in this work. One-way surface chemical reactions are considered in the first model. The rates of backward reactions 
are assumed to be negligible. This one-way reaction model is obtained by modifying the carbon surface ablation model 
developed by Park[7]. In the original Park’s model, there were not enough surface reactions to account for surface 
catalysis and thus the surface heating was under-predicted. Additionally, its nitridation reaction does not appear to 
agree with available data. In our previous work, the nitrogen recombination reaction was proposed, and the rate of 
carbon sublimation was adjusted to be consistent with measurements from the arc -jet tests. The carbon surface is 
assumed to be fully catalytic for nitrogen recombination reaction. 

The second model is more complex compared with the first one. It contains two-way surface reactions. This model 
is derived from the approach developed by Zhluktov and Abe. [8] Zhluktov’s model includes surface catalysis, 
oxidation, and sublimation at the carbon surface. The reactions considered in the original model by Zhluktov and Abe 
were found to be inadequate to predict both surface heating and surface recession because some key reactions were 
not considered. Two additional reactions are included in this work: one is for the oxidation reaction between molecular 
oxygen and the carbon surface, and the other is the recombination of nitrogen atoms. The forward and backward 
reaction rates computed in our study use the general formulation developed by Marschall,[13] in which the forward 
reaction rates (Kf) were written in 5 different forms. They are Eley-Rideal, Sublimation, Adsorption, Arrhenius, and 
Langmuir-Hinschelwood. The proposed gas/surface reactions and their forward reaction rate constants are listed in 
Tables 1 and 2: 


GSI Model I: 





Kf 

Reaction Type 

Source 

1 0 + C(b) CO 

V 4 J 

x 1.0 


Eley-Rideal 

Park(99)[7] 
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2 0 2 + 2C(b) -4 2CO 


x0.5 


Eley-Rideal 


Park(99)[7] 


3 2N -> N 2 


^ xl.O 
4 


Eley-Rideal 


CG(13)[5] 


4 3C(b) -4 C 3 


x (8.27 x 10 14 ) x e 


9.084 wo*) Sublimation 


CG(13)[5] 


5 C 3 ^ 3C(b) 


x0.25 

4 


Eley-Rideal 


Park(99)[7] 



Reaction Type Source 


1 O + E(s) O(s) 


x 1.0 

40, 


Adsorption 


ZA(99)[8] 


2 20(s) ~ 0 2 + 2E(s) 


3.583 x 10 xT xe 


3 . 079 x 1 o 4 ! Arrhenius 


ZA(99)[8] 


3 0 2 + E(s) O + O(s) 


Eley-Rideal 


ZA(99)[8] 


4 0 2 + E(s) + C(b) «-» CO + O(s) 


— x 0.5 x e~ T 


1 . 44 M 0 5 1 Eley-Rideal 


5 O(s) + C(b) <-> CO + E(s) 


2.080x10' xT xe 


-3.999x10*1 Arrhenius 


ZA(99)[8] 


6 O + O(s) + C(b) C0 2 + E(s) . i 2 . 0 xi o 5 

I v n ) T 


Eley-Rideal 


ZA(99)[8] 
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7 20(s) + C(b) <-> C0 2 + 2E(s) 


3.583 x 10 15 xT xe 


3.999S<10 ) 

T 


Arrhenius 


ZA(99)[8] 


8 C + E(s) E(s) + C(b) 


9 C 2 + 2E(s) <-» 2E(s) + 2C(b) 


10 C 3 + 3E(s) 3E(s) + 3C(b) 


11 N + E(s) ~ N(s) 


12 N 2 + E(s) ~ N + N(s) 


13 N(s) + N(s) N 2 + 2E(s) 


f — \ 
v c 

1 4 ®J 


x 1.0 


f tt A 




x 1.0 


fV| 

1 


xO.25 


f ~ \ 

v N 

1 4 ®J 


xl.O 


f rr A 


v 4 ® s y 


x e 


7.659M0* 


f 


323 


n_ 6.0221x10 


Eley-Rideal ZA(99)[8] 


Eley-Rideal ZA(99)[8] 


Eley-Rideal ZA(99)[8] 


Adsorption ZA(99)[8] 


Eley-Rideal ZA(99)[8] 


Langmuir-Hinschelwood * 


In the above proposed reactions, C(b) represents the solid carbon species, E(s) is the numbers of empty surface 
sites, and O(s) and N(s) are sites occupied by oxygen and nitrogen, respectively. In computing the forward reaction 

rate constants,® , is the active site density, and v k is equal to -yj8RT w / JlM k . Backward reaction rate constant can 

be calculated from chemical equilibrium constant and forward reaction rate constant for each two-way reaction. 

The DPLR code was enhanced by Mac Lean based on Marschall’s formulation to solve the surface species mass 
balance equation for the finite-rate gas surface reactions. The detail of MacLean’s implementation can be found in 
Ref. 9. All flow simulations presented in this work are performed using this version of DPLR code (DPLR/4.03.1). 

There are twenty two gas phase chemical species used in this study for the simulation of PICA and air-Argon 
interactions. None of the ion species are considered. The chemical species are C 0 2 , CO, N 2 , 0 2 , NO, C 2 , C3, CN, H 2 , 
HCN, C, N, O, H, CH, CH 2 , C 2 H, C 2 H 2 , C4, C3H, C4H, and Ar. The gas phase chemical reactions and their rates 
adopted in the DPLR simulations are listed in Table 3. 
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K f = AT n e~ Ta/T 

A, cc/mol/s 

n 

T n ,K 

Source 

1.* 

C 0 2 + M <-> CO + 0 + M 

6.9xl0 21 

-1.5 

63275 

PH(94)[14] 


M = C, N, O 

1.4xl0 22 

-1.5 

63275 

PH(94)[14] 

2* 

CO + M<->C + 0 + M 

2.3xl0 20 

-1.0 

129000 

PH(94)[14] 


M = C, N, O 

3.4xl0 20 

-1.0 

129000 

PH(94)[14] 

3.* 

n 2 + m~n + n + m 

7.0xl0 21 

-1.6 

113200 

PH(94)[14] 


M = C, N, 0 

3.0xl0 22 

-1.6 

113200 

PH(94)[14] 

4.* 

O2 + M 0 + 0 + M 

2.0xl0 21 

-1.5 

59750 

PH(94)[14] 


M = C, N, O 

l.OxlO 22 

-1.5 

59750 

PH(94)[14] 

5.* 

NO + M<->N + 0 + M 

5.0xl0 15 

0 

75500 

PH(94)[14] 


M = C, N, O 

l.lxlO 17 

0 

75500 

PH(94)[14] 

6.* 

c 2 + m<-*c + c + m 

3.7xl0 14 

0 

69000 

PH(94)[14] 

7. 

C3 + M<->C 2 + C + M 

6.3xl0 16 

-0.5 

101200 

BL(92)[15] 

8 . 

CN + M<->C + N + M 

2.5xl0 14 

0 

71000 

PH(94)[14] 

9. 

h 2 + m<->h + h + m 

2.2xl0 14 

0 

48300 

SB(93)[ 16] 

10. 

NO + O ~ O2 + N 

8.4xl0 12 

0 

19450 

PH(94)[14] 

11. 

N 2 + O <-* NO + N 

6.4xl0 17 

-1.0 

38370 

PH(94)[14] 

12. 

CO + 0 <-» O2 + c 

3.9xl0 13 

-0.18 

69200 

PH(94)[14] 

13. 

CO2 + 0 <-» O2 + CO 

2.1xl0 13 

0 

27800 

PH(94)[14] 

14. 

CO + C C2 + 0 

2.0xl0 17 

-1.0 

58000 

PH(94)[14] 

15. 

CO + N CN + 0 

l.OxlO 14 

0 

38600 

PH(94)[14] 

16. 

N 2 + C <-> CN + N 

l.lxlO 14 

-0.11 

23200 

PH(94)[14] 

17. 

CN + 0 NO + C 

1.6xl0 13 

0.10 

14600 

PH(94)[14] 

18. 

CN + C C 2 + N 

5.0xl0 13 

0 

13000 

PH(94)[14] 

19. 

HCN + H <-> CN + H 2 

8.0xl0 n 

0.0167 

85537 

He(67)[17] 

20. 

CH + M<->C + H + M 

1.9xl0 14 

0 

33700 

DH(92)[18] 

21. 

CH 2 + M « C + H 2 + M 

1.3xl0 14 

0 

29700 

DH(92)[18] 

22. 

CH 2 + M <-» CH + H + M 

4.0xl0 15 

0 

41800 

DH(92)[18] 

23. 

c 2 h + m«->c 2 + h + m 

1.74xl0 35 

-5.16 

57400 

KR(97)[19] 

24. 

C2H2 + M <-> C 2 H + H + M 

6.96xl0 39 

-6.06 

67130 

KR(97)[19] 

25. 

C 2 + C 2 *■* C3 + c 

3.2xl0 14 

0 

0 

KR(97)[19] 

26. 

C 2 + H 2 «-* C 2 H + H 

6.6xl0 13 

0 

4030 

KR(97)[19] 

27. 

CH 2 + C <-> CH + CH 

1.62xl0 12 

0 

23600 

M(67)-NIST[20] 

28. 

ch 2 + ch 2 « c 2 h 2 + h 2 

1.58xl0 15 

0 

6010 

Ba(95)-NIST[20] 

29. 

CH 2 + C C 2 H + H 

5.0xl0 13 

0 

0 

DH(92)[18] 

30. 

CH 2 + C 2 H CH + C2H2 

1.81xl0 13 

0 

0 

TH(86)[21] 

31. 

CH + CH ^ C 2 H + H 

1.5xl0 14 

0 

0 

DH(92)[18] 

32. 

CH + C 2 H C2H2 + C 

l.OxlO 14 

0 

0 

DH(92)[18] 

33. 

CH 2 + CH 2 C2H2 + H + H 

2.0xl0 14 

0 

5530 

Ba(95)-NIST[20] 

34. 

C2 + C2 + M C4 + M 

l.OxlO 27 

-3.0 

0 

KS(92)[22] 

35. 

C + CH C 2 + H 

5.0xl0 14 

0 

0 

KR(97)[19] 

36. 

C + c 2 h <-> c 3 + h 

7.0xl0 14 

0 

0 

KR(97)[19] 

37. 

C 2 + CH C 3 + H 

2.0xl0 14 

0 

0 

KR(97)[19] 

38. 

C 2 + C 2 H C 4 + H 

1.2xl0 14 

0 

0 

KR(97)[19] 

39. 

CH + H <-» C + H 2 

l.lxlO 14 

0 

0 

KK(93)[23] 

40. 

CH 2 + H ~ CH + H 2 

l.lxlO 14 

0 

0 

KK(93)[23] 

41. 

C 2 H + C 2 H «-» C2H2 + C 2 

l.OxlO 13 

0 

0 

KS(92)[22] 

42. 

C + C2H2 <-» C3H + H 

4.0xl0 14 

0 

6000 

KK(93)[23] 

43. 

C 2 + C 2 Hz *-* C 4 H + H 

l.OxlO 14 

0 

0 

KS(92)[22] 

44. 

C2 + C4H C2H + C4 

1.2xl0 14 

0 

0 

KS(92)[22] 
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45. 

C 2 H + C 2 H C 4 H + H 

l.OxlO 14 

0 

0 

KS(92)[22] 

46. 

C 3 H + H <-> C 3 + H 2 

l.OxlO 14 

0 

0 

KK(93)[23] 

47. 

C 4 H + H <-> C 4 + H 2 

2.0xl0 13 

0 

0 

KS(92)[22] 

48. 

H + C 4 + M <-> C 4 H + M 

1.74xl0 37 

-5.5 

0 

KS(92)[22] 

49. 

CH + CH C 2 + H + H 

l.OxlO 14 

0 

0 

KK(93)[23] 

50. 

C 2 H + H 2 ~ C 2 H 2 + H 

7.4xl0 14 

0 

3350 

KR(97)[19] 


* two temperature reaction rates 
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III. Surface Thermal Chemistry Coupling Simulation 


The finite-rate gas-surface interaction models proposed in this work are implemented by means of the surface 
thermal chemistry coupling between DPLR and TITAN. The interface between a TPS material and its surrounding 
flow field can be defined by solving species mass conservation and energy balance equations. Species mass 
conservation at the surface of TPS material is written as: [24] 


-pDjVCt + pv w C, = Wj + m C, 


'g i>8 


(i) 


The first term on the left-hand side is mass transfer through diffusion, and the second term is mass transfer due to 
convection. On the right-hand side are the source terms due to gas-surface interaction and pyrolysis gas injection. 
Based on global mass balance at the surface, the following equation for the total mass blowing rate is expressed as: 


pv w = m c + m g 


( 2 ) 


The total convective heat flux to the surface for the flow field that includes a two -temperature model is given as: 


Vconv = ~ k t VT t ~ k v VT v 


JhiP D y C i- 


( 3 ) 


Energy conservation equation at the surface is written as 

Rconv 3" ^w) "f" ^g(j^g ^w) tt w q rw (J£ W (T W fa, ) tfcond 0 (4) 

The first term in Equation (4) is the total convective heat flux, the second and third terms represent the heat of ablation, 
the fourth and fifth terms are radiation absorption and emission, respectively, and the final term is the rate of heat 
conduction into the TPS material. 

To obtain the solutions for Equations (1) to (4) requires computations of thermal and species diffusion rates in the 
flow field at the surface as well as thermal diffusion and pyrolysis gas injection rates of TPS material at the surface. 
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In our DPLR/TITAN coupled simulation, the governing equations for fluid and solid are solved separately, and 
Equations (1) to (4) have to be solved either in the flow code or in the material response code. Thus, communication 
between two codes needs to be established for exchanging information on surface thermal chemistry and shape change. 
For a general finite -rate surface boundary condition, the species mass conservation equation is solved with flow-field 
governing equations and the total energy balance equation is solved with solid material governing equations. The 
detail implementation of surface species mass conservation equation (Eq. 1) solved with DPLR was discussed in the 
work of MacLean. Surface thermal chemistry and shape change have to be shared between the flow solver (DPLR) 
and the material solver (TITAN) using either explicit or implicit approach. The schematic diagram in Figure 1 depicts 
how this coupled simulation is performed. The pyrolysis gas is assumed to be at chemical equilibrium before being 
injected into the adjacent air. This assumption was proved to be reasonable for conditions studied in this paper. [25] 
The computational grids for both flow field and TPS material are reconstructed based on the shape change predicted 
by the material response code, and the free stream conditions are also updated according to the flight trajectory. For 
arc -jet simulations, the free stream conditions remain unchanged during the exposure time. 



Figure 1. Surface thermochemistry and shape change coupling methodology. 
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IV. Results 


The ablator simulations presented in this section include two parts. In the first part, we discuss the computation 
of stagnation models. One group of stagnation models is with relatively low arc stream total enthalpy, between 5 and 
12 MJ/Kg, for which tests were conducted in the Aerodynamic Heating Facility at NASA Ames Research Center [26]. 
The other group of stagnation models with relatively high stream total enthalpy from 15 to 30 MJ/kg, for which tests 
were conducted in the Interactive Heating Facility at NASA Ames Research Center [27]. In the second part, the 
simulations of SPRITE test model, which is a 55 -deg half -angle spherical cone with stream total enthalpy ranging 
from 5 to 26 MJ/kg, are presented. The computations in the first part focus on the prediction of total stagnation 
recession, and those performed in the second part focus on the analysis of both surface recession and shape change. 
The GSI models discussed early in this paper are examined and the recession predictions are compared with the 
measurement to understand the accuracy of these proposed GSI models. The coupled DPLR/TITAN simulation system 
developed in our previous work is used to perform all the computations presented in this section. 

Part 1 : 

The purpose of performing stagnation model simulation is to examine the proposed carbon/air interaction 
models, which include a one-way reaction model (Model I) and a two-way reaction model (Model II). The arc -jet tests 
selected for computations in Part 1 have stagnation heat fluxes from 45.2 to 1100 w/cm 2 , which is sufficient to cover 
the rate controlled and diffusion controlled carbon surface oxidation, as well as carbon sublimation. Five of the test 
conditions, with relatively low total enthalpy (below 12 MJ/kg), selected from the Mars Science Laboratory (MSL) 
test series have stagnation heat fluxes of 45.2, 72.8, 114, 133, and 236 W/cm 2 . The other seven high enthalpy test 
conditions (above 15 MJ/kg ) with heat fluxes of 107, 169, 246, 395, 552, 694, and 1 102 W/cm 2 are obtained from the 
Orion Thermal Protection Systems Advanced Development Project (TPS ADP) test series. The arc -jet test stream 
conditions of these selected cases are listed in Tables 1 and 2 for the low enthalpy MSL and the high enthalpy TPS 
ADP, respectively. Thus, the oxygen molecules may not be fully dissociated in the low enthalpy MSL arc -jet stream, 
and are fully dissociated in the high enthalpy TPS ADP cases. Figure 2 shows the shape of stagnation model used in 
this part of simulation. The nose radius equals to the model diameter, and the sides are cylindrical. This is the so called 
“iso-q” geometry (r n = 5.08 cm, and r c /r n = 1/16). The test material is PICA, and the model holder is made of LI-2200. 
There is an air gap with the size of 1.5 mm in-between the back face of PICA and the front side of LI -2200. 
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Two-block grid system is used for both DPLR (99 x 159 and 23 x 159) and TITAN (120 x 60 and 40 x 80) 


simulations. The grid adaption option is turned on to perform grid alignment with the bow shock wave. The boundary 
condition index used in DPLR allows the user to specify surface species mass injection rate, temperature, and finite- 
rate surface reactions. The implicit coupling approach between DPLR and TITAN is used to predict PICA surface 
recession. [6] 


Table 1: Low enthalpy MSL arc -jet cases selected for analysis. 


Case number 

Stagnation point heat flux, 
W/cm 2 (cold wall) 

Stagnation point pressure, 
kPa 

Exposure time, 
s 

Total enthalpy 
MJ/kg 

1 

45.2 

4.88 

240 

5.5 

2 

72.8 

8.85 

120 

6.4 

3 

114 

20.9 

80 

5.9 

4 

133 

31.6 

80 

6.1 

5 

236 

25.0 

35 

11.1 


Table 2: High enthalpy Orion TPS ADP arc-jet cases selected for analysis. 


Case number 

Stagnation point heat flux, 
W/cm 2 (cold wall) 

Stagnation point pressure, 
kPa 

Exposure time, 
s 

Total enthalpy 
MJ/kg 

1 

107 

2.3 

55 

15.2 

2 

169 

5.0 

60 

17.0 

3 

246 

8.5 

42 

19.3 

4 

395 

17.2 

34 

21.4 

5 

552 

27.3 

30 

23.3 

6 

694 

31.0 

27 

29.2 

7 

1102 

84.4 

10 

25.6 



Figure 2. Geometry and material map of model used in the arc -jet tests. 
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The comparison of predicted stagnation point recession with arc -jet data for the low enthalpy (MSL) group is 
shown in Fig. 3a. The circles in red represent the post arc -jet test measurement. Three different models of surface 
chemistry are used for recession predictions. The first one is the chemical equilibrium surface model (diamonds) [27], 
the second one is Model I with one-way finite -rate GSI (triangles), and the third one is Model II with two-way finite- 
rate GSI (squares). The chemical equilibrium surface model is established based on the pre-generated B’ tables for 
PICA exposed to air/ Argon mixture. The chemical equilibrium model over -predicts the recession by about 35% for 
the heat flux condition of 45 W/cm 2 , and by about 25% for the heat flux condition of 73 W/cm 2 . This result is expected 
because carbon surface oxidation for such low surface heating conditions is not likely to be in chemical equilibrium 
and needs to be determined by the finite -rate computation. For the heat flux conditions of 114 and 133 W/cm 2 , the 
recession predictions using chemical equilibrium model are about 20% higher than the arc -jet data, and for the 236 
W/cm 2 case, chemical equilibrium recession prediction is about 5% higher than the data. Generally speaking, the 
recession predictions using Model I are slightly higher than those using Model II. Model II appears to agree the best 
with arc -jet data among three different surface chemistry models, especially at the lowest heat flux condition, 45.2 
W/cm 2 . Chemical equilibrium surface model is commonly used for the recession prediction of many arc -jet test 
simulations because it is the simplest and least expensive approach, but also it produces reasonable recession 
predictions for most cases. In the lowest heat flux condition, where the backward reaction rates can no longer be 
ignored, a two-way reaction model thus has to be considered. In Fig. 4a, predicted temperature and pressure 
distributions along stagnation streamline for condition with heat flux of 45.2 W/cm 2 are presented. The stagnation 
pressure is near 5 KPa. The maximum translational temperature is around 4500 K, and that of vibrational temperature 
is 4000 K. Figure 4b shows the distributions of selected species mass fraction (N 2 , O 2 , NO, CO, and O) along 
stagnation streamline. The black lines are the predictions using the gas/surface interaction Model II. The red lines are 
those predictions using a chemical equilibrium surface assumption, which is achieved by increasing the reaction rates 
six orders of magnitude. The surface recession computed using this approach is the same as that computed using the 
B’ tables[27]. The near surface CO concentration for chemical equilibrium model is higher than that for finite-rate 
model. Consequently, the near surface O 2 concentration for chemical equilibrium model is lower than that for finite- 
rate model. All of the atomic oxygen is consumed at the carbon surface. 
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Figure 3a Comparison of computed stagnation point 
recession with arc -jet data for low enthalpy MSL cases. 
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Figure 3b Comparison of computed stagnation point 
recession with arc -jet data for high enthalpy Orion 
TPS ADP cases. 


Similar comparison of computed stagnation point recession with arc -jet data for the high enthalpy group (CEV 
ADP) is presented in Fig. 3b. The chemical equilibrium model has been presented in the work of Milos[27], and it 
was shown that the predictions using chemical equilibrium model agreed reasonably well with the data for these 
selected test conditions. In the first six conditions, the carbon surface oxidation is controlled by the oxygen diffusion 
rates, and in the seventh condition, the carbon recession is dominated by the process of sublimation. For the diffusion 
controlled conditions, both Model I (triangles) and Model II (squares) give almost identical predictions, and agree 
with arc -jet data (circles). For the condition with the highest heat flux, 1 100 W/cm 2 , Model I prediction is around 5% 
below data, and Model II prediction is 20% lower than the data. In model II, the backward reaction rates of carbon 
sublimation are computed from chemical equilibrium constants and forward reaction rates. Even if we tune up the 
forward rate constants of sublimation reactions, the predicted net total recession still remains about the same because 
the net computed reaction rates are limited by chemical equilibrium constants. The reason for Model II under- 
predicting the sublimation rate by about 20% needs to be investigated further. There is no oxygen involved in the 
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sublimation process. The sublimation rate also depends on the concentrations of gas phase C 3 , C 2 , and C at the surface. 
Thus, the accuracy of predicting carbon sublimation rate is directly tied to the accuracy of both gas phase and gas- 
surface chemical reaction models. The central focus of current study is not the gas-gas chemical reactions. In our 
future work, the study of gas phase chemical reaction models will be addressed to further understand and improve the 
prediction accuracy of sublimation rate described in Model II. 



Figure 4. a: Pressure and temperature distributions 
along stagnation streamline for stagnation heat flux of 
45.2 W/cm 2 . 



Figure 4.b: Comparison of mass concentration 
distributions along stagnation stream for stagnation 
heat flux of 42.5 W/cm 2 . 


Part 2 : 

The SPRITE arc -jet test model studied in this part was used for the investigation of a conformable ablative TPS 
technology targeting missions requiring peak heat-flux around 250 W/cm 2 . The model is a 20.32-cm (8-inch) diameter, 
55-deg half-angle sphere-cone with a 5.08-cm (2-inch) spherical nose radius. The heat shield consists of 3 segments 
of conformable PICA and 1 segment of standard MSL PICA. All segments have a thickness of 1.27-cm (0.5-inch). 
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The back shell is entirely LI-2200. This work focuses on the segment with MSL PICA material, because its material 
properties are well studied and thus is a good candidate material for GSI model validation. 

The cross section of SPRITE test article is shown in Fig. 5. a. Simulation is performed by assuming that the test 
model remains as an axi-symmetric system during the entire exposure time. The nose cap (in orange) is made of PICA, 
and the segment (in yellow) is also made of rigid PICA. There is a gap between these two elements. The rest of region 
is made of LI-2200. The back face of PICA segments and the gaps between segments were filled with RTV. However, 
RTV is not considered in this simulation for simplicity. The initial axi-symmetric computational grids used for 
simulation are presented in Fig. 5.b. A one -block grid system is used for DPLR (133 x 99) and a two-block grid system 
is used for TITAN (102 x 153 and 31 x 123). Both grids are moving with time to reflect the surface recession. 
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Figure 5.b: The initial axi-symmetric computational grids 
for DPLR and TITAN used in the coupled simulation. 


Figure 5. a: The cross section of a SPRITE test model. 


Four arc -jet test conditions for SPRITE model are considered for analysis. Their conditions are listed in Table 3. 
The highest stagnation-point cold-wall heat flux is 700 W/cm 2 , and the lowest one is 90 W/cm 2 . The heat flux at mid- 
conical section is about 50 to 60% of that at the stagnation point. The stream total center line enthalpy is in the range 
from 5.4 MJ/kg to 26.0 MJ/kg. For total enthalpy at 5.4 MJ/kg, the oxygen molecules in the arc flow are not fully 
dissociated. Additionally, the surface temperature for this lowest enthalpy condition, especially in the conical section 
(flank), should remain low and the carbon oxidation is controlled by the reaction rates. 
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Table 1: SPRITE arc -jet cases selected for analysis. 


Case number 

Stagnation point heat flux, 
W/cm 2 (cold wall) 

Stagnation point pressure, 
kPa 

Exposure time, 
s 

Est. total enthalpy 
MJ/kg 

1 

700 

31.5 

30 

26.0 

2 

340 

17.9 

60 

17.0 

3 

290 

9.1 

80 

18.0 

4 

90 

11.7 

100 

5.4 


The PICA surface recession and shape change predicted by the coupled DPLR-TITAN simulation system are 
presented in Figs. 6. The gray lines represent the model initial shape before its exposure to the arc stream. The red 
lines are the post-test measurement of model surface shape. The predicted model surface at the end of exposure is in 
blue for Model I (one-way reaction model) and is in black for Model II (two-way reaction model). One -dimensional 
point calculation with the assumption of chemical equilibrium surface is also performed at the stagnation point and 
the mid-point of conic section (flank). The chemical equilibrium computations are calculated using FIAT and the 
results are shown as the black dots in Figs. 6. For Conditions 1 to 3, the blue lines (Model I) are on top of the black 
lines (Model II). Thus, the black lines are hardly to be seen. This indicates both models give almost identical 
predictions and agree very well with the post-test measurement. The chemical equilibrium predictions are generally 
slightly highly than those predicted by the finite-rate models, and are still in good agreement with measurements for 
Conditions 1 to 3. The carbon surface for Conditions 1 to 3 is in the diffusion controlled oxidation regime. All three 
surface chemistry models give consistent predictions with the measurement. Based on our predictions, it can be 
concluded that the complicated finite-rate simulation does not seem to be necessary, and a simple chemical equilibrium 
model should be adequate to predict the SPRITE model shape change with reasonable accuracy if carbon oxidation is 
limited by the diffusion rate. On the other hand, if the surface heating is low, such as that for Conditions 4, chemical 
equilibrium surface model over-predicts the surface recession by about 20% at the stagnation point, and 35% at flank, 
as shown in Fig. 6.d. The shape change and recession predicted by using Model II (black line) agree well with the 
measurement, except the areas covered by glass fence around the outer edges of gaps. These glass fences were 
generated from the product of decomposed RTV, as the surface temperature is not high enough to reach the melting 
point of glass (-1600 °C). Glass accumulates on the carbon surface near gaps to form these fences. The formation of 
glass fences may affect the accuracy of the axi-symmetric computation model. However, it should have little impact 
in the nose cap region because the subsonic flow is confined to a small region near stagnation point. Thus, the 
stagnation heating should not be affected by the presence of downstream glass fences. Besides, if the majority of solid 
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glass accumulation occurred later in the exposure time, the computation error associated with the presence of glass 
fences should be even lower. These fences are not seen in the other conditions with higher surface heat fluxes because 
their surface temperatures are higher than the glass melting point. For Condition 4, the recession prediction on the 
surface of the nose cap using Model I is also in good agreement with data, but Model I over-predicts the recession by 
about 20% at flank. It indicates that under such conditions the backward reaction rates become significant. This finding 
is consistent with the results obtained for stagnation models discussed in Part I. 
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Figure 6. a: Comparison of surface recession and shape 
change between computations and measurement for 
condition 1. 


Figure 6.b: Comparison of surface recession and shape 
change between computations and measurement for 
condition 2. 




Figure 6.c: Comparison of surface recession and shape 
change between computations and measurement for 
condition 3. 


Figure 6.d: Comparison of surface recession and shape 
change between computations and measurement for 
condition 4. 
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Part 3: 


The uncertainty on the arc stream total enthalpy is estimated to be within 15-20%, surface heat flux is 10%, and 
pressure is 5%. Additionally, the uncertainty of the surface recession measurement is around 10%. Consequently, we 
estimate that the error on the surface recession prediction should not exceed 15%, if the surface chemistry is correctly 
modeled. 

In Fig. 7a, surface pressure is plotted against surface temperature for all arc -jet conditions studied in this work. 
Diamonds represent the low enthalpy MSL conditions, squares represent the high enthalpy CEV ADP conditions, and 
circles represent the SPRITE stagnation point and mid-frustum conditions. For conditions in which the predictions 
using chemical equilibrium surface chemistry agree with measurements within ±15%, their symbols are colored in 
blue. If the agreement between chemical equilibrium surface prediction and measurement is greater than ±15% and 
that between finite-rate surface prediction and data is less than ±15%, their symbols are in red. Thus, for those 
conditions in red, the chemical equilibrium surface assumption is unlikely to be true and the finite -rate model has to 
be employed to correctly predict carbon recession. In Fig. 7a, it can be seen that, at a given surface pressure, the 
conditions at higher temperature are in blue and those at lower temperature are in red. This is because the reaction 
rates increase with increasing temperature. Figure 7b shows a similar chart for surface pressure vs. stream total 
enthalpy. All the symbols in red are at the low enthalpy side, and those in blue are with relatively high enthalpy. It 
indicates that surface pressure is not a key factor to determine the surface chemical state as compared with enthalpy 
and temperature. Figure 7c shows the chart for total enthalpy vs. surface temperature. At a given surface temperature, 
the conditions with higher enthalpy are at chemical equilibrium, but those with lower enthalpy are not. This suggests 
that total enthalpy is the most important factor to determine the chemical state of surface for arc stream conditions 
studied in this paper. Figure 7d presents the chart for boundary edge oxygen molecule concentration vs. total enthalpy. 
It can be found that the O 2 concentrations for symbols in red are significantly high than those in blue. O is far more 
reactive than O 2 when they interact with a carbon surface. This can explain why a carbon surface exposed to high 
enthalpy flow is more likely to reach a chemical equilibrium state. 
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Figure 7a: Arc -jet test model conditions, surface 
pressure vs. surface temperature. 


Figure 7b: Arc=jet test model conditions, surface 
pressure vs. stream total enthalpy. 
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Figure 7c: Arc -jet test model conditions, stream total 
enthalpy vs. surface temperature. 


Figure 7d: Arc -jet test model conditions, boundary 
layer edge oxygen molecule mass fraction vs. stream 
total enthalpy. 
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V. Conclusions 


Two sets of finite-rate gas-surface interaction models for carbon based ablator have been proposed. The predicted 
PICA surface recession using the proposed models with surface heat fluxes ranging from 40 to 1100 W/cm 2 were 
compared with arc -jet data. The predictions using the model with two-way surface reactions agree well with data if 
carbon recession is primarily an oxidation process. However, the proposed two-way reaction model, presented here, 
under-predicts the carbon sublimation rate. The predictions using the proposed one-way reaction model are also in 
good agreement with data except the relatively low heat flux condition (below 50 W/cm 2 ), in which the backward 
reaction rate becomes significant. Generally speaking, the chemical equilibrium surface model tends to over-predict 
the carbon surface recession, if carbon oxidation is rate controlled, and a finite-rate model has to be applied to 
accurately compute the carbon recession. The simulation of SPRITE arc -jet model reconfirms that the proposed two- 
way reaction model can accurately predict the recession and shape change of PICA material at low surface heat flux 
condition, where the chemical equilibrium and one-way reaction models are normally over-predicting. It was also 
found that the carbon ablation process is more likely to reach a chemical equilibrium state if the surface is exposed to 
a high enthalpy flow with fully dissociated oxygen, as compared to a low enthalpy flow with partially dissociated 
oxygen. 
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